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Abstract. We classify four-state spin models with interactions along the edges 
according to their behavior under a specific group of symmetry transformations. 
This analysis uses the measure of complexity of the action of the symmetries, in 
the spirit of the study of discrete dynamical systems on the space of parameters 
of the models, and aims at uncovering solvable ones. We find that the action 
of these symmetries has low complexity (polynomial growth, zero entropy). We 
obtain natural parametrizations of various models, among which an unexpected 
elliptic parametrization of the four-state chiral Potts model, which we use to 
localize possible integrability conditions associated with high genus curves. 
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1. Introduction 

Our purpose is to describe the structure of the set of spin models with interactions 
along the edges, and propose a classification of the matrix of Boltzmann weights of 
such models. 

The main tool we use is the action of a group of transformations which was 
exhibited in previous works concerning integrable models ||l|, ||, H) , and which has its 
origin in the so called inversion relations Q. They are generated by basic involutions, 
realized as birational transformations of the parameters of the models (we will call 
this group of transformations I-symmetries) . 

We analyze the behaviour of the realizations of I-symmetries, especially the 
iterates of specific infinite order elements of the group. The guiding principle is 
that a chaotic behaviour of the realization of I-symmetries, signaled by a measure of 
their complexity (entropy) , is preventing integrability of the corresponding statistical 
mechanical model. 

We concentrate on the four state models because of their richness, but the method 
applies to all values of the number q of states. 

The paper is organized as follows: we describe the models we have in mind 
and the action of I-symmetries. We then describe the classification of the models 
according both to the dimension of their parameter space and to their behavior under 
I-symmetries. We recover, in the course of this classification, a number of known 
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two-dimensional models which have been studied in the literature, among which the 
standard scalar Potts model the Ashkin- Teller model the Kashiwara-Miwa 
model 0, and the four-state chiral Potts model ||^. One section is devoted to the 
analysis of definite models. Among other results, we give an unexpected^ elliptic 
parametrization of the chiral Potts model. We give a constructive way to write 
explicitly the algebraic locus where the group of I-symmetries degenerates into a finite 
group, leaving room to solvability. We finally comment on the specificity of four-state 
models as compared to three and five or more states. 

2. Spin edge models, I-symmetries, and admissible patterns 

The models are defined on a D-dimensional regular lattice (square, triangle, 
honeycomb, cubic, . . . ), each vertex bearing a spin a which can take q values. A model 
is defined by giving the matrix of Boltzmann weights of local edge configurations: the 
state of an edge is determined by the values of the spin at its ends. One arranges the 
weights in a matrix W whose entry Wij is the weight of the edge configuration with 
one end in state i and the other end in state j. Notice that edges are oriented, and 
W is not necessarily symmetric. One may furthermore distinguish different types of 
bonds, like, for example, vertical and horizontal ones on a square lattice, and the most 
general non isotropic model is then associated to a number of q x q matrices Wi, , one 
for each type of bond. 

Since the entries of the matrices W are Boltzmann weights, W is defined up to 
an overall multiplicative factor. The most general matrix is then a collection of q^ 
homogeneous entries, defined up to a factor. 

Simple transformations may be defined, which act on the matrices W, forming 
the group of I-symmetries. This group is generated by simple involutive generators, 
namely the matrix inverse /, and the element by element inverse J, both taken up to 
a overall multipHcative factor. They may be written 

/ : Wki~^ Aki (1) 

J : Wki~~^ l/Wki (2) 

where Aki is the cofactor of Wki in W. To have a well defined action of J, we assume 
that entries are non vanishing. 

The two generators are non commuting involutions. They generate an infinite 
group isomorphic to the dihedral group Z2 x Z. Its infinite part is made of the 
iterations oi ip — I ■ J and its inverse — J ■ I. 

We will build families of models by successive reductions of the number of 
parameters, obtained by imposing equality relations on the entries of the Boltzmann 
matrix W. 

To illustrate this, consider the chiral Potts model: its Boltzmann weight matrix 
is a cyclic matrix W, where instead of keeping all entries independent, one imposes 
the following equality conditions: 

Wu,v = Wu+l,v + l (3) 

I Especially_Li|i view of the occurrence of higher genus curves in the solutions of the star-triangle 
relations M, nd] 
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getting for g = 4 a pattern of the form: 

/Wo 

W = 

V 
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Wl 


W2 




W3 


Wo 


Wl 


Wl 


W2 


ws 


Wo 



(4) 



instead of the general 4x4 matrix. The number of parameters has gone down from 
16 to 4 because we have imposed 12 equaUty relations. 

The specificity of this type of relations is two-fold. Any relation of the form 
Wa,b — Wa'fi' for some pairs of indices {a, 6}, {a', 6'} defines a pattern of matrix 
which is automatically left stable by the generator J of the group of I-symmetries. 

Some pattern defined in this way will also be left stable by the matrix inversion 
/, hut not all of them. We will say that a pattern is admissible if its form is left stable 
by matrix inversion ||l|. 

Giving a pattern is equivalent to giving a partition of the entries of W . The 
number of possible patterns, that is to say the number P{q^) of partitions of objects, 
is given by 



q s-1 



,s=i k=a 



k) 



k\ (s - 1 - fc)! 



(5) 



It grows extremely rapidly with the number q of spin states. It is 21147 for g = 3, 
10480142147 ~ 10^° for g = 4, and 4638590332229998592 ~ 4.6 • lO^^ for g = 5. 

Our first step is to determine which of the P{q^) patterns are admissible. For 
q = A this was done by direct inspection of the P(16) ~ 10^° patterns. The outcome 
is a list of 166 patterns. We have not considered vanishing conditions on the entries, 
and we restrict ourselves to invertible matrices ||, so that both / and J are defined 
and invertible. 

There are trivial redundancies in this list because any of the q\ permutations of 
the q values of the spin, i.e. the same permutation acting on simultaneously on the 
rows and columns of W will take an admissible pattern into an admissible pattern: it 
acts by similarity on W , and commutes with I-symmetries. This reduces the list to 



only 42 different admissible patterns, given in Appendix A . 

We see that there are no admissible patterns with 9, 11, 12, 13 , 14, or 15 free 
(homogeneous) parameters. 

A matrix of the form @ of the previous section belongs to the similarity class of 
pattern #17, which contains three elements (see later): 



/Wo Wl W2 W3\ 



W3 
W2 



Wo 

W3 
W2 



Wl 
Wo 
Ws 



W2 
Wl 



/Wo Wl W2 W3\ 



,Wi W2 W3 Wo Ws W2 Wl Wo/ \W2 W3 Wi Wo. 

These patterns have four homogeneous independent parameters [wq, wi,W2, w^]. Since 
they are admissible, there is an action of the homogeneous matrix inversion / on 
[wo,wi,W2,W3], which reads, when written for the first of the three representatives 
(cyclic matrix): 



W2 
Wl 



Wo 
W3 
W2 



W3 
Wo 
Wl 



Wl 
W2 



/Wo Wl W2 waX 



Wl 
W-i 
W2 



Wo 
W2 



W3 

Wo 



W2 
Wl 



Wowl 



W2w1 



Wo — > Wq — 2 W0W1W3 
Neither have we considered equality up to sign between entries | ^ 



w\w2 
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Wi > — + 2 WQW1W2 + W1W3 — W3W0 ~ ""^2^3 (6) 

W2 — » ^2 — 2 i(;ii(;2W3 — W2W0 + ^'^oWi + wIwq 
wz — > — + 2 W0W2W3 + w^wf — w\w\ — WqWi 
The action of J reads, if we use the homogeneity of the entries: 

[W0,W1, W2, W3] »[W1W3W2, WoW3W2, W0W1W3, Wi W2 Wq] 

3. Complexity analysis 

Since I-symmetries are essentially the iterates of f ^ I ■ J, their action looks like a 
discrete time dynamical system defined on the parameter space of the model. The 
measure of complexity we use is given by the rate of growth of the degrees of the 
iterates of (p. Indeed once we have reduced the action of / and J to the pattern, 
the action of is a polynomial transformation on the independent parameters. All 
expressions have a definite degree d. The n-th iterate naively has degree d" , but since 
we work with homogeneous coordinates, we should factor out any common factor, so 
that the degree may drop to a lower value d„. 

A measure of the complexity of the map is given by the entropy 

e= lim -log(d„). (7) 

n — ^oo Ji 

If e = the growth of c?„ is polynomial in n and d„ behaves like 

dn = an'' {l + 0{l/n)) (8) 

and one may then define a secondary complexity index associated to the map, 
the integer power k. 

The most compact way to encode the sequence of degrees d„ is to write its 
generating function: 

n—oo 

gis) = ^ rf„s". (9) 

11=0 

This function was conjectured to always be rational with integer coefficients, i.e. 
be of the form: 

9is) = §4 (10) 

Q[s) 

with P and Q some polynomials of degree p and q with integer coefficients. It is 
sufficient in practice to calculate the p + q first terms of the sequence to infer this 
generating function[|. Any further degree calculation serves as a verification. 

If the roots of the denominator of the generating function are all of modulus one, 
then the entropy vanishes, and we may consider the secondary index k. 

To illustrate the method, we may look at pattern #17 as in the previous section. 
The degree of ip is 9. The first iterations of (p yield the sequence of degrees: 

[1, 9, 33, 73, 129, 201, 289, 393, 513, 723, . . .] (11) 
The beginning of the sequence (up to dg) may be fitted with the generating function 

{l + 3sf 

5(5) = 7- ^ (12) 

(1 - s) 

II The catch is that the values of p and q are not known in advance. 
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which is compatible with the above-mentioned conjecture. The subsequent terms are 
then be predicted correctly, and gn is found. We see that the growth of the degree is 
polynomial (e — 0) and the secondary index is k = 2. 

For each pattern the different representatives yield the same transformations and 
consequently the same generating function. 

If the entropy e vanishes, the secondary index may be read from the generating 
function g: if g{s) has a pole of order 7 at s = 1, then k + 1 = 7. 

We have written the realizations of (/? = /• J for the admissible patterns, and 
calculated the generating function for all of them. The result is that for all 42 
admissible patterns the growth of the degree is polynomial, i.e. the entropy vanis hes, 
and the inde x k takes values 0, 1, 2 and 3 (see the following sections, Appendix A and 



Appendix B] ). 



The complexity can also be estimated using the following arithmetical procedure: 
start with a matrix with random integer entries, and iterate ip, dividing all the 
coefficients by their greatest common divisor. The complexity can be estimated from 
the growth of a typical entry. In practice the number of digits of the entries grows like 
A", with e = log A. Since the entries of the matrix grow fast, one has to use special 
representation of the integers allowing to manipulate arbitrary large values |l4|. In 
practice one performs the iteration with several initial random integer matrices. The 
average l{k) of the logarithm of the entries is recorded as a function of the number of 
iterations k. The entropy e or the secondary index k are easily deduced from l{k). 



4. Structure of the set of admissible patterns 

There exists a partial order relation on the set of patterns, induced by the partial 
order on the partitions of the entries. Indeed a partition a = {ai, a2, ■ ■ ■ , a^} may be 
finer than a partition /? = {f3i,(32, ■ • ■ if all parts constituting a are subsets of 

the constituents of /3: 

a^P, if Vfc = l.../i, 3/ s.t. afe C A (13) 

If a /? will say that /? is a descendant of a. We obtain a descendant by adding 
further equality relations on the entries, i.e. merging parts together. 

As an illustration, pattern #3 i.e. the standard scalar Potts model is a descendant 
of pattern #17 i.e. the chiral Potts model. All patterns are descendant of pattern 
#42 which has the finest partition. 

The list of immediate descendants of the various patterns can be given with the 
notation [ Pattern number, { list of direct descendants}], where {} indicates an empty 
list. 

[4, {2}], [5,{1}], [6,{}], [7, {1,3}], [8, {3}], [9, {2}], [10, {4, 5}], 

[11, {4, 6}], [12, {4}], [13, {6}], [14, {5}], [15, {7}], [16, {5, 7}], 

[17, {6, 7}], [18, {9}], [19, {8}], [20, {5, 6}], [21, {12}], [22, {8}], 

[23, {4, 7}], [24, {9}], [25, {15}], [26, {16, 17, 20, 25}], 

[27, {13, 14, 20}], [28, {19, 22}], [29, {12, 15, 23}], [30, {10, 14, 16, 23}], 

[31, {11, 13, 17, 23}], [32, {18, 24}], [33, {10, 11, 20, 21}], 

[34, {8, 9, 23}], [35, {21, 25, 29}], [36, {23}], 

[37, {26, 27, 30, 31, 33, 35}] [38, {18, 22, 29, 34, 36}], [39, {24, 30, 34}], 
[40, {26, 30, 36}], [41, {19, 31, 34, 36}], [42, {28, 32, 37, 38, 39, 40, 41}] 
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Figure 1: Diagram of descent relations. 
5. Extended permutation classes 

If two patterns Hi and 112 verify 112 = L • Hi • i?, where both L and R arc permutation 
matrices, the realizations of (pi = Ii ■ Ji (resp. (p2 = I2 ■ J2) we get from Hi (resp. 
112) are again of the same degree, even if L and R are not inverse of each other. We 
even have ip\ = ip2- This explains why in the table of generating functions, non similar 
patterns may share the same generating function. In the figures, we have used framing 
boxes to indicate these equivalences, listed in the following table: 



# param. 

2 

3 

4 



6 
7 



Classes 

ri = {1,2,3} 

r2 = {4,5,6,7} 

r3 = {8,9}, r4 = {io,i7} 

r5 = {ll,16}, r6 = {12,13,14,15} 

r7 = {18,22}, r8 = {19,24} 

r9 = {20,23}, rio = {21,25} 

Tn = {26,30,31,33}, ri2 = {27,29}, 

ri4 = {34,36}, ri5 = {35} 



ri3 = {28,32} 
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8 


Tie 


= {37} 


10 


ri7 


= {38}, ri8 = {39,40,41} 


16 


Tig 


= {42} 



Remark: Lp^ commutes with matrix transposition t as well. A a consequence, t 
acts on the classes T. It actually acts trivially on these classes. 



6. Gauge equivalence 

It is important to take into account the usual gauge symmetries of the Boltzmann 
weights. Their action may be represented by a similarity transformation: 







W 



W' g-^ - W ■ g with g 






Vo 



S2 






S3 





(14) 



Si, 



They commute with the generator / of I-symmetries, but not with J. However, if two 
matrices verify jl^), then their images by J are automatically gauge related by 

/1/si 



9 



\ 



1/S2 
1/S3 






This means that J has a well defined induced action on the gauge equivalence classes, 
although it does not commute with the gauge transformations. Transformations as 
( p^ yield the following equivalences: 

Te = {12, 13, 14, 15} = {4, 5, 6, 7} 

r7 = {18,22} _-^r3 = {8,9} 

ri2 = {27,29} — .r9 = {20,23} 
ri3 = {28,32} — .r8 = {19,24} 



Through these equivalences we have reductions by one unit of the number of 
parameters. Moreover for all these reductions, there exists an /-invariant which can 
be changed by the action of gauge transformations, and the reduction amounts to fixing 
the value of this invariant. Take for example pattern #12. One of its representatives 
has Boltzmann matrix: 



Wi2 = 



/Wo 


Wl 


W2 


W2 \ 


Wl 


Wo 


W2 


W2 


W3 


W3 


Wl 


Wo 


\W3 


W3 


Wq 


Wl / 



The quantity: 

2 \W3 W2 
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is /-invariant. It can be brought to the value A = 1 by a gauge transformation of the 
form (|lj) with si — S2 , S3 = S4, and u>2 s§ = 82- This makes W[2 — ■ W12 ■ g 
symmetric and yields pattern #4, with Boltzmann matrix: 



W[2 ^Wa^ 



/Wo 


Wi 


W2 


W2 \ 


Wi 


Wo 


W2 


W2 


W2 


W2 


Wi 


Wq 


\w2 


W2 


Wo 


Wl / 



We may picture the gauge equivalences by the following diagram: 




Figure 2: Gauge reductions 



The usual models which have been studied in the literature are present in the top 
of this Hst (they have a low number of parameters) . The standard scalar Potts model 
is pattern #3. The symmetric Ashkin- Teller model is pattern #7. The Ashkin- Teller 
(alias Kashiwara-Miwa |0| with = 4) model is pattern #16. The chiral Potts model 
is the pattern #17. 

They all share the property that I-symmetries have rational invariants, i.e. the 
orbits of the realization of I-symmetries are confined to algebraic varieties. Since we 
know how crucial is the role played by these varieties in the solvability of the models, 
it is natural to look for the structure of invariants of the admissible patterns. 
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7. Invariants of the action of I-symmetries 

Since all elements in an extended class yield the same birational realizations, up to a 
permutation of coordinates, invariants may associated to the classes T . 

The descendance of a given pattern being obtained by adding relations among 
the entries of W, descendants inherit induced invariants (rational or not). Some 
invariants will disappear or coalesce in the process. What might also happen is that 
a transcendent invariant degenerates to a rational one. 

The value of k is related to the existence of conserved algebraic varieties under 
the action of I-symmetries. If there exist invariant algebraic curves, k is necessarily 
0,1 or 2 H, If the curves are rational, k = or 1. The values of k are shown on 
the following diagram. 




K=2 



34) 



38 



35 



37 

K — J 



K=3 



39 40 



m 



Figure 3: Values of the index k 

All the values of n are corroborated by the procedure presented in paragraph 3. 
A number of iterations of the order of thirty and an averaging over a number of the 
order of ten initial matrices is enough to get a good precision on the value of k. 

There exists an unbounded algorithm to find rational invariants ||l2|. There is 
however no simple way to discover the non-rational ones, and it is an open problem to 
find all /-invariants. When the number of parameters is small, a graphical approach 
may be extremely useful (see later with pattern #25). In the next section we will 
examine in detail the structure of typical patterns. 
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8. Analysis of specific patterns 

8.1. Pattern #3; Standard scalar Potts model 
The Boltzmann matrix is: 



In term of the parameter w = Wi/wq, the action of ip = I ■ J reduces to a homographic 
transformation: I : w ^ —2 — w, .J : w ^ ^/w, (p : w —2 — 1/w. In terms of the 
variable t = 1/(1 + w), ip is the translation t — > t — 1. 



/Wo 


Wl 


Wl 


Wl\ 


Wi 


Wo 


Wl 


Wl 


Wl 


Wl 


Wo 


Wl 


\wi 


Wl 


Wl 


Wq/ 



8.2. Pattern #7; Symmetric Ashkin-Teller model 



The Boltzmann matrix 


is: 












(Wo 


Wl 


W2 


W2\ 




Wr = 


Wl 
W2 


Wo 
W2 


W2 
Wo 


W2 
Wl 






\w2 


W2 


Wl 


Wo J 




and there is one algebraic invariant 






WqWi 

^7 = , 

W2 [Wl 


2 

- W2 
- Wo) 






The curves Ay = a have a rational parametrization: 


Wo = t, 


Wl 


1 


— at 


W2 = 1 


t 


— a 


8.3. Pattern #8 












The Boltzmann matrix 


is: 












/Wo 


Wl 


Wl 


Wl \ 




W8 = 


Wl 
Wl 


W2 
W3 


W3 
W2 


W3 
W3 






\wi 


W3 


W3 


W2 / 




The quantity 












ZAg — 


W2 Wl 


'^ + 2wi'^W3 - 


2 Wo W3 W2 






W2 Wi'^ - 


- Wo Ws'^ 



(15) 



^ (16) 

is invariant by J and changes sign under the action of /. The action of ip exchanges 
the surfaces S+ (Ag^'' = a) and E_ (Ag^' = —a). Since the condition Ag^'' = a may 
be solved rationally in Wo, we may write the action of (p'^ on S_|_, with coordinates 
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wi,W2,W3. In terms of the inhomogeneous variables x = wi/w^ and y = ^3/(^2+^3), 
it reads: 

, {b + by-y) 



(2 + y) {by~2y + b- 

y — >y' = y + b 

with 6 = 2a/(l + a), and the quantity 



1) 



.(2) 



X with 



$(y) ^ r 



(17) 
(18) 

(19) 
(20) 



b J \ b J' \ b 
is invariant by the action of ip'^ on E+. This example, where k = 1, possesses a 
mixture of algebraic and non-algebraic invariants, which can be evaluated exactly. 
The orbits of (p are confined to non- algebraic curves, and the n-th iterate may be 
written explicitly: 

Hyo) 



yn = yo + nb, 



Xo 



(21) 



Hyn) 

The orbits accumulate to the point (xoo,yoo) — (Ag (xo, yo), 00), which is a point on 
the line ^2 + ^3 = 0. 

Pattern number 8 leads to pattern number i if — wi and W2 — wo, in which 
case the above two invariants take definite values: Ag^'' = 00 and Ag^' = 1. 



8.4- Pattern #16; Ashkin-Teller model 

This model coincides with the iV 4 Kashiwara-Miwa model 
matrix is: 

/Wo Wi W2 W3\ 
■Wi Wo W3 W2 



16 



W2 W3 Wo Wi 



The Boltzmann 



(22) 



> W3 W2 Wl Wo . 

There are two independent algebraic /-invariants 



^16 



Wl Wz — Wo W2 



A 



(2) _ Wl W3 - Wo W2 



(23) 



Wo Wl — W3 W2 Wo W3 — Wl W2 

Fixing the values of these two invariants to oi and 02 determines an elliptic curve in 
the space of parameters. The equation of this curve is 



aix^ — aiy^ — a2X — a2aiy + a2y^x 



with — 

W3 



W2 _ Wo 

W3 ^' W3 



a2aiyx = 
X + aiy 
y + aix 



(24) 



8.5. Pattern #17; Chiral Potts model 

The chiral Potts model is equivalent to pattern #17. Its Boltzmann matrix was already 
given above (equation (ji])). There exist two algebraically independent invariants. 

(1) W3 W2 - Wo Wl 



A 



Potts 



Wl W2 — Wo W3 



^(2) ^ 
^Potts 



{wo W2 — Wl^) {wo W2 — W3^) 
{wo W3 - Wl W2)^ 



(25) 
(26) 
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The condition Ap^^^^ = a may be solved rationally in wo- Setting ^pl^^ = together 
with X = wi /w2 and y — /w2 , yields the equation of an eUiptic curve to which the 
orbit is confined: 

{x — y)^ {x + y)^ b — {—yx^a + ax + x^ — y) {—y + ax + y^x — y^a^ — (27) 

We may illustrate the efficiency of the graphical method in this case, since the set 
of the iterates of ip is dense in a curve^. Figure 4 shows a typical orbit in the three 
dimensional space of parameters. The orbit is projected on a coordinate plane. The 
aspect of the picture is very stable under changes of the starting point of the iteration. 




Figure 4: A typical orbit oi tp — I ■ J for the chiral Potts model 



The transformations / and J may be restricted to the quadrics a 
may be written on x and y, for fixed a. 

"■^ — 2 y'^xa + 2ya + x'^ya^ — a^x — x 



y 



1 



y 

J(i . X 



y = 



y^x 



ay"^ + a'^yx 
2 x^ay — y 



xy — a"^ + ax"^ 
ya^ + 2ax + x^c? 



1 



ay"- 



a^yx 



xy — a'^ 

x' = l/x, y — >y' = l/y 
1 Contrary to what happens in the example of pattern #8 above. 



ax^ 



est. They 
(28) 

(29) 
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The transformations la and Ja leave the curves ( |27| ) globally invariant whatever a and 
b are. 

The modular invariant of the elliptic curve (|2^) is given by 

52 = 1 + 256 b'^ - 512 - 16 b^a* + 288 b^a^ - 16 6^ + 16 a^b 

- 32 a^fo + 16 a6 - 4 + - 4 + 6 a'^ 
33 = (l - 32 6^ + 32 afe - 2 + a^) (l + 16 6^ - 16 afo - 2 + a") 

X (1 - 8 6^ + 8 a6 - 2 + a'') /(3V3) 

The modular invariant (|o|) can be written 
. _ (i-A^+M2)3 

where 

(l-a)2(l + a)2 ^ ^ 

is the square of the modulus of the elliptic functions parametrizing the curve. Notice 
that A4 is ambiguously defined by condition (^ij) , and may be replaced by any of the 
six values {M, 1 - M, l/M, 1 - 1/M, 1/(1 - M),M/iM - 1)}. 

There are equivalent symmetric biquadratic forms of the curve (|2^ 

_ 2 {J,\jy^ + Jy^J,^ + Jjj,^)pq + 4 Jjjy^j,^ ip + q) + 

~^ (^Jx Jy ~^ Jy Jx Jz ) ^ Jx Jy Jz i^Jx ~^ Jy ~^ Jz ) — 

and 

{Jx - Jy) [p'q^ + 1) - {Jx + Jy) (p' + g") + 4 J,pg = 

with 

J 2 _ J 2 



Jn,^2a, Jy^4:b-2a, J,^a'^ + 1, and M 



,y — -lU ^ U,, ,^2 — " I J-, ailU J^l — 2 _ T 2 

Changing 7W to 1 — tW or 1 /M amounts to permuting Jx, Jy, and Jz ■ 

Generic orbits of (p — I ■ J are infinite, but we know how important are the 

degenerate cases where these orbits are finite: for two-dimensional lattice models, this 

is where non trivial star-triangle integrability takes place |^, 17 1. For given values of 

a and b, the action of (/5 is a shift on the curve (p7|). 

It is possible to write down the Weierstrass form of the curve 

- a X - + ^0 (32) 

with0 

" ' /3=^7r^ (33) 



48(6-a)8' 288 (6-a)i2 

and give an explicit coordinate transformation from the original variables (x, y) to 
[X,Y]. We obtained the transformation through Maple's implementation of the van 
Hoeij algorithm |[l^ , after setting b = a + l/u^. 

+ The normalizations could be changed, and the powers of l/(f)— a) could be absorbed in the definition 
of X and Y. 
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The action of (p, when written on ( |3^ is the addition^ of the point 11+ with 
coordinates [X^^^ , Yy^^]: 

1 + 16 62 + 8 a6 - 2 a2 + 12 6 + 12 0^6 + a* 
- - 1/12 -4 

(b~a) 



b(l + ar h+4:b-2a + a^) 

Y —1/2 — ^ - '- 

(b-a) 

Although the change of coordinates depends on u, the [X, Y] coordinates of the point 
n+ depend only on and may be written rationally in terms of a and b. 

The inverse map is the addition of the opposite point n_ = [Xt^^^j—Yt^^^]. 
Notice that both 11+ and n_ are obtained from the double point x = —l/au,y ~ 00 
of ( p^ ) in the desingularisation process which leads to the Weierstrass form (|3^). 

The sum of n± with itself, denoted 2 • Il± , is given by 

[-1/12 l±i^^±Hi^l^iJ^^^±^, ±1/2 + (34) 
(6 — a) (6 — a) 

and the multiples n ■ Il± are then easily obtained by recurrence. Their [X, Y] 
coordinates depend only on and may be re-expressed rationally in terms of a and 
b. 

The point 2 • 11+ corresponds to the point 

(l + a^jau (2 + au^) a^u 
1 — M^a'^ — 1 — u a — 



of the original curve ( P7|) . 

The condition = ic? is just: 

X{n-U+) = X{m-n-) and y(n • n+) = -r(m • n_) 

As an example, the condition that </? is of order 3 is: 

V3^{l + af b-a^ ^0 (35) 

which is an algebraic condition on the entries of the Boltzmann matrix, and can be 
checked directly. 

The condition that (p is of order 4 is the vanishing of F(2 • 11+) given in (|^ i.e. 

V4 + (2 6 -a) a 

The condition 2 6 — a = is precisely the condition appearing in |^ |l^, where 
solutions of the star-triangle relations related to higher genus curves were exhibited. 
The condition ^3 = is not known to yield similar integrability conditions and was not 
found to yield any transition in the statistics for the Random Matrix Theory approach 
(see II). 



8.6. Pattern #25 

This example is a good illustration of the mixed use of algorithmic research of algebraic 
invariants, graphical method, and the knowledge of the entropy and index k. The 

To add two points on the cubic @, draw the straight line through the two points. Compute the 
third point of intersection of the line with the curve, and take its symmetric under Y —Y. 
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Boltzmann matrix is: 



W- 



25 



/Wo 



Wl 
Wo 
W3 



W2 
Wq 



W3^ 
W2 

J 



,W4, W2 W4 Wq, 

The number of inhomogeneous parameters is four, 
algebraic invariants. 



. (1) _ Wq Wl - W3 W2 
WqWs — Wl W2 



,(2) 
^25 



W4 W3 — Wo 



'■) ( 



Wl W4 — W2 



(36) 



We found two independent 



(37) 



(38) 



(W4 W3 — U'2^) {Wi W4 — Wo"^) 

Drawing orbits of the iterates of if leads, with a proper choice of starting point, to 
pictures like Figure 5. Contrary to what happens with pattern #17 (Figure 4), the 
shape of the orbit is rather unstable under changes of the starting point. 




Figure 5: A typical orbit of /-symmetries for pattern #25 

Such a picture indicates that there exists another independent invariant. Would that 
invariant be algebraic, the stable curves would be elliptic, and we know that the index 
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K would then be less than 3. Since k = 3 we know this additional invariant has to be 
non-algebraic. 

Equation of the surface Si : Ajs^ = a may be solved for wi. We get an expression 
of / and J on Si. Changing variables to 

[x,y,Z,t\ = [wl,WoW2,WoW3„WzW4\ 

we obtain: 

Ia{[x, y, z, t]) ^ [x (x + ay - at - tf (x + y) {x + ay) , 

{x + y) {—yx + xta + xt — y'^a^ {x + ay) {x + ay — at — t) , 

— z {x + y) (x^ + xta'^ — xt ~ y'^a'^') {x + ay ~ at ~ t) , 

{x + ay) {x-y)t {x^ + xta^ - xt - y^a^)] (39) 

with remaining algebraic invariant 

X (t ~ y) (yt — yx — y^a + axt) ^ ^ 

A^5 = :f '- (40) 

(x-y) t{x + y) 

Notice that in equation (|39|), the third coordinate z is just multiplied by a factor 
depending on the other coordinates, and does not appear anywhere else in the induced 
transformation la nor in the invariant A25. The curves A25 = b are eUiptic curves 
in the variables {x,y,t), extending to cylinders in [x,y,z,t]. The curve appearing on 
Figure 5 is drawn on such a cylinder. 

We may examine a particular value of b where the curve A|g = b degenerates 
to a rational curve, but there exists a transcendent invariant for la- The rational 
parametrization of the cylinder is 

,s'^t(l + sa) st(l + sa) , , , 

[x,y,z,t] = [ 4 ^,z,t 41 

s + a s + a 

The induced map (/ • J)^ may be written 

S^,^S, ^^^^ (42) 

with T{S) = ^— 5— and 

(1 - S^)(l - S) 



s - g-i , ^^ a+l + v/(l-a)(l + 3a) 
S= , q=-l/2 ^- . 



s-q 

If we introduce the infinite product 



k — OG 



then the 2nth iterate of (i.j) reads 

S2n = g'" ^0, Z2n - ^0 ' (44) 

i.e. the transcendent quantity A — z ■ Il{S) is invariant by (/ • J)^. The convergence 
of these Eulerian products is ensured when a £ (—1/3,1). Such values of a would 
produce orbits with accumulation points. If a is outside this interval, then q has unit 
modulus and the orbits are similar to the one shown in Figure 5. The quantity n(S') 
is replaced by the standard analytic prolongation. 
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9. Other values of q 

We have investigated the structure of the set of admissible patterns for other values 
of q. 

The value g = 4 is an important threshold. It is known to be one for the standard 
scalar Potts model in statistical mechanics and for the chromatic polynomials 
and coloring problems. What we show is that it is the maximal value ensuring the 
vanishing of entropy, whatever admissible pattern we choose. 

In the case of a cychc matrix of size q (chiral Potts model), with entries 
[wq, . . . , Wq-i] there is a simple relation between matrix inverse and element by element 
inverse. Define the linear transformation C (discrete Fourier transform, or Kramers- 
Wannier duality) 

9-1 

C : Wk — > '"^Wr, 

where uj is the g-th root of unity. We then have: 

i = c-^ - J -C. 

It is possible to analyze the singularity structure of the transformations / and J, and 
prove that the generating function for the degree of the iterates of Lp for the chiral 
Potts model with q states is ||l^ 

(l-5)(l-s(g2-4g + 2) + 52)- 

As soon as g > 5, the entropy becomes strictly positive for this family of admissible 
patterns. 

For all sizes q, there exists a pattern depending on q{q + l)/2 homogeneous 
parameters: the symmetric matrix. 

An open question is to decide whether or not there exist patterns depending on 
r parameters, with q{q + l)/2 < r < q^. For g = 3 and q = 4 there is a no man's land 
between q{q + l)/2 parameters (symmetric matrix) and q^ parameters (generic q x q 
matrix) . 

Another open question is to know if there is anything between cyclic and 
symmetric matrix when q is prime. 

The value q = 2 is trivial and yields a finite group of symmetries. When (7 = 3, 
the same analysis can be performed completely (see below). When q = 5, the number 
of patterns is too large to be thoroughly examined, and we have restricted ourselves 
to patterns with only three homogeneous parameters. 

Of course there exist admissible patterns with larger values of q still yielding 
vanishing entropy. They may have for example a block structure, the blocks being 
themselves admissible patterns of smaller sizes pOl El], p3, p3|. 



9.1. Three state models 

The analysis can be made along the same lines as for the four-state models: 
enumeration, calculation of the entropy, calculations of invariants, etc. It is 
summarized here. 

The number of partitions of nine objects is 21147 and the number of admissible 
patterns is only 9, when equivalences are taken into account. The list of patterns is 
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given in [Appendix C . For patterns #1, #2, #3, #4 the group of I-symmetries is finite. 
Patterns #3 and #4 are related by an extended permutation similarity as defined in 
section (||). Patterns #5 and #6 yield the same generating function for the degrees 
of the iterates of Lp 

^ (l + 2s)(l + s)^ 

(1 + s + .s2) (1 - sy 

and the trajectories of (p are confined to algebraic curves. Patterns #7, #8, #9 yield 
the generating function 

1 + 6 s + 9 s2 + 2 + 6 .s^ + 3 s'^ - 6 - 3 , , 

97 = .98 = .99 = -; "2 ■ (46) 

(1 + .S + S2) (1 - S f 

Patterns #7 and #8 (pattern #8 is the 3x3 symmetric matrix) are also related by an 
extended p ermutation si milarity. A list of algebraic invariants is given for the various 
patterns in Appendix D . 

For pattern ^^8 (symmetric 3x3 matrices), the Boltzmann matrix is 

/ Wo Wi W2 ' 



Wl W3 Wi 



(47) 



V W2 Wi W5 , 

If we normalize the entries with the condition wq = \, the action of p? reduces to the 
following homothetic transformation: 



with 



[l,Wl,W2,W3, W4, Ws] > [1, Ai Wl, A2 W2, A3 W3, A4 W4, A5 W5] 

(W5 - W2^) W3 W4 (-Wl W2 + W4) 



A2 = 



(-W3 W5 + W/^)W2 (-W2 W3 + Wl W4) 
(W3 - Wi^) W5 W4 (-Wl W2 + W4) 



( — W3 W5 + W42) (— Wl W5 + W4 W2) Wl 

A3 = A?, A4 = Ai A2, As = A2 

where Ai and A2 are left invariant by the action of ^p. 

The appearance of third root of unity in ^ has to do with the specificity of the 
third power of Lp. 

There is a parametrization of the general three-state model (pattern #9) which 
simplifies the action of I-symmetries. The matrix of Boltzmann weights reads 

W = 

If one introduces the quantities 

2 a 2 2 

a = Wq , /) = W4 , 7 = w's 

d — w\ W3, e = W2 We, = wr 

Wl We W5 

A = — , ^ = — , iy= — 

W3 W2 W7 

The array [a, (3, 7, 5, e, Q transforms under pr' by multiplication of its entries by factors 
[ka, ki3,kj,ks,ke, kc]. One can check that all the ratios ka/kp, . . . are invariants of ipfl . 
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In a similar way, the three quantities [A, /i, v] transform under (^^ by multiplication by 
factors [kx, fc^, ki,] which are all invariants of ip^ . 

This yields a system of algebraic invariants of rank four of ip^ , for example: 



W2 W3 wr ' Wi'^we W2 ' wi W2 W3 wq ' W2 we wr 
9.2. Five state models and beyond 

The number of partitions of 25 objects is over 4.10^*. To enumerate the configurations 
is beyond reach. We have restricted our analysis in the following way: if one 
looks for partitions with only three parts, that is to say five-state models with 
only three possible colors (three homogeneous parameters), then there are "only" 
141197991025 ~ 141.10^ patterns. We have performed their full enumeration, and 
found the subset of admissible ones. The interesting result is that there exits a unique 
admissible pattern, up to permutations of rows and columns: the symmetric Potts 
model, with symmetric cyclic 5x5 matrix of Boltzmann weights 



/ Wo Wi W2 W2 Wi \ 
Wi Wq Wi W2 W2 



w = 



1^2 Wi Wq Wi W2 
W2 W2 Wi Wo Wi 
\wi W2 W2 Wi Wq/ 

There are actually 72 representatives, all related by independent permutations of rows 
and columns, as in section (||). 

The trajectories of the iterations of ip are confined to algebraic curves. The 
entropy vanishes and the secondary index k has value 2 | |l5| , [l^ . The corresponding 
statistical mechanical model is known to be integrable, via star-triangle relation, 
when these curves degenerate to rational curves and the symmetry group becomes 
finite |, 0,0, 

One should keep in mind that admissible patterns of size 5x5 easily yield 
transformations with exponential growth of the degree. For example the cyclic matrix 
(chiral Potts model, see eq. (Eq)) yields a map with non vanishing entropy 



1 /7 + 3\/5, ,^ , 

e = log( ^ ), (48) 

and the corresponding statistical mechanical model is generically not star-triangle 
integrable. 

We have used the arithmetical method described in section (^) to evaluate the 
entropy for both the symmetric pattern (symmetric 5x5 Boltzmann matrix) and 
the general five-state model (5x5 Boltzmann matrix with no relations between the 
entries)|| . We found the same value for the entropy as for the cyclic matrix, given by 
equation (|4^). 

In contrast, if one restricts the 5x5 Boltzmann matrix to be both cyclic 
and symmetric, the entropy vanishes, and the transformation ip has one algebraic 
invariant [ p5[ . 

One may wonder if the entropy calculated for the cyclic matrix, the symmetric 
matrix, and the most general qx q matrix still identify for q > 6. One may also try to 
evaluate the entropy for the cyclic-symmetric pattern, and provide a closed expression 
valid for all values of q, and similar to (0). 

(t For a generic 5x5 matrix the coefficients got as large as lo^^'""^" after six iterations. 
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10. Conclusion and perspectives 

We have given a classification of four-state spin models according to the properties of 
the realization of I-symmetries. 

This preparatory work provides the optimal parametrizations of the models. 

If one is looking for integrability via solutions of star-triangle equations, one 
should stick to cases where the parameters live on invariant algebraic varieties. Those 
may be present for generic values of the parameters, or appear by a collapse of 
the realization of the infinite group to a finite group, or by the degeneration of 
transcendental invariants to algebraic ones. 

The algebraic approach presented here may be complemented with random matrix 
theory, as in ||l9|, to locate solvability of the lattice models. 

One may have to examine in detail the properties of both the matrix product 
and the element by element product, and not just the stability by the corresponding 
inverses, when mixed with matrix transposition. One should notice that the 
Boltzmann matrices of the models which are known to be integrable via the star- 
triangle relation, i.e Fi, with the notations of section ||, are (bi)-stochastic 
matrices (sum of elements on lines and columns are constant) . 

A detailed analysis of the star-triangle equations will then proceed on a case by 
case basis, which is beyond the scope of this paper. 

Appendix A. List of the 42 admissible patterns for q — ^ 

The sixteen entries of each 4x4 matrix are presented on one line in the order 
[Wii,Wi2,Wi3,Wi4,W2i,W22, ■ ■ ■ ,W4:3,W4i]. The choice of the representative is 
arbitrary, and should be understood modulo permutations of Hues and columns. The 
list is arranged with increasing number of independent homogeneous parameters, and 
this number is indicated after the list of entries. The last item of each line is the 



number 


of patterns in the similarity class. 














#1: 


[wo 


Wo, Wo, wi 


Wo 


Wo, Wl 


Wo 


Wo, Wl , Wo 


Wo 


wi,wo, Wo 


Wo] 


2 , 


3 


#2: 


[wo 


Wo 


Wl , Wo 


Wo 


wi,wo. Wo 


Wi,Wo, Wo 


Wo 


Wo, Wo, Wo 


Wl] 


2 , 


6 


#3: 


[wo 


Wl 


Wl , Wl 


Wl 


wo,wi, Wl 


wi,wi, Wo 


Wl 


Wl,Wl, Wl 


Wo] 


2 , 


1 


#4: 


[wo 


Wl 


W2,W2 


Wl 


Wo, W2 


W2 


W2, W2 


Wl 


Wo 


W2,W2, Wo 


Wl] 


3 , 


3 


#5: 


[wo 


Wl 


W2,Wo 


Wl 


Wo, Wo 


W2 


W2, Wo 


Wo 


Wl 


Wo,W2, Wl 


Wo] 


3 , 


3 


#6: 


[wo 


Wo 


Wi,W2 


Wo 


Wo, W2 


Wl 


W2, Wl 


Wo 


Wo 


Wl,W2, Wo 


Wo] 


3 , 


3 


#7: 


[wo 


Wl 


W2,W2 


Wl 


Wo, W2 


W2 


W2, W2 


Wo 


Wl 


W2,W2, Wl 


Wo] 


3 , 


3 


#8: 


[wo 


Wl 


Wi,Wi 


Wl 


W2, W3 


W3 


Wi,W3 


W2 


Ws 


Wi,Ws, W3 


W2] 


4 , 


4 


#9: 


[wo 


Wl,Wl, Wl 


Wl 


W2, W3 


W2 


Wi,W3 


W2 


W2 


Wi,W2, W2 


Ws] 


4 , 


12 


#10 


■ [wo 


Wl 


W2,W3 


Wl 


Wo, W3 


W2 


W2, W3 


Wl 


Wo 


W3, W2,Wo 


Wl] 


4 , 


3 


#11 


■ [wo 


Wl 


W2,W3 


Wl 


Wo, W3 


W2 


W3, W2 


Wl 


Wo 


W2,Ws, Wo 


Wl] 


4 , 


3 


#12 


■ [wo 


Wl,Wl, W2 


W3 


W2,Wo, W3 


W3, Wo 


W2 


Ws 


W2,Wi, Wl 


Wo] 


4 , 


3 


#13 


■■ [wo 


Wl 


Wo,W2 


W2 


W3 , Wl 


W3 


Wo, W2 


Wo 


Wl 


Wi,Ws, W2 


Ws] 


4 , 


3 


#14 


■■ [wo 


Wl 


W2,Wo 


Wl 


W3, W3 


W2 


W2, W3 


W3 


Wl 


Wo,W2, Wl 


Wo] 


4 , 


3 


#15 


■■ [wo 


Wl 


Wi,W2 


Ws 


Wo, W2 


W3 


Ws, W2 


Wo 


Ws 


W2,Wi, Wl 


Wo] 


4 , 


3 


#16 


■ [wo 


Wl 


W2,W3 


Wl 


Wo, W3 


W2 


W2, W3 


Wo 


Wl 


W3, W2,Wi 


Wo] 


4 , 


1 
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#17 




Wo 




Wi 


Wn 


Wq 


Wey 


11) Q 


Wo 


Wn 




Wo 


UJS 


Wl 


\ 


#18 


f . u; / 


Wo 


Wn 


W<i 


W/ 
4 


W-Q 


W-Q 


Wo 
^z 


Wi 


Wn . 


Wn 


Wn 


Wl 


Wn 


Wo] 


#19 




Wo 


W-Q 


Wl 


W/ 
4 


Wi 


Wl 


W<i 


Wi 


Wn . 


Wo 


Wo 


Wl 


W-Q 


Wn] 


#20 




We) 


Wn 


Wq 


Wn 


Wn 


W 1 


W 1 


Wn 


ll)n 


^0 






11)1 


tlJn] 


#21 


\wn Wi 


id 


Wi 


Wq 


Wp 


W J 


Wn 


Wo 
^z 


Wl 


Wn 


Wl 


W 1 

4 


Wn 


Wq 


Wo] 


#22 


[if n . Wi 


Wi 


Wi 


Wo 


Wq 


W 1 
4 


W / 


Wo 


W 1 
4 


Wq . 


W / 
4 


Wo 
^Z 


W / 

^4 


W 1 


Wq] 


#23 


\wn . / 


Wo 


W^ 


Wl 


Wn 


Wo 
^z 


W0 
^z 


Wo 
^z 


Wo 
^z 


W-Q . 


W) 
4 


Wo 


Wo 


W / 
4 


Wq] 
"-'0 J 


#24 


[w)n Wi 


Wi 


Wi 


Wl 


11)0 


Wq 


W 1 
4 


Wi 


Wq 


H) 1 
^4 ' 


Wo 


Wl 


W4 


Wo 
^z 


IIJq] 
ujs \ 


#25 




Wo 


W-Q 


W/ 
4 


Wn 


W/ 

4 


Wo 


Wo 
^z 


W<i 


Wn . 


Wl 


W/ 
4 


Wo 


W/ 
4 


Wn] 


#26 


[ "-^ (/ 7 ^ i 


We} 


Wq 




Wo 


Ws 


W5 


Ws 


Ws 


Wo, 


W4 


W3 


W2 


Wl 


Wn] 


#27 


N/M 7/;, 


Wn 




W3 


W4 


W5 


W4 


Wq 


W2 


Wo, 


Wl 


Ws 


W4 


Wq 


W J ] 


#28 


\inn 11)1 


W 1 


Wi 


W2 


Ws 


W4 


W5 


W2 


Ws 


Ws, 


W4 


W2 


W4 


Wc: 


^Dq] 
UJS \ 


#29 


\wn . W 1 

[ f "■'1 


Wo 


Wi 


Ws 


W4 


Ws 


Ws 


W2 


Wl 


Wo, 


Wl 


Ws 


W5 


W-Q 


wA 
UJ4 J 


#30 


\wn . W 1 


Wo 


W-Q 


Wl 


W4 


Ws 


Ws 


W2 


Ws 


Wo, 


Wl 


Ws 


Ws 


Wl 


W/] 
UJ4 J 


#31 


Wi 


Wo 


Wq 


Ws 


W4 


Wl 


Ws 


W2 


Ws 


Wo, 


Wl 


Wl 


W5 


Wq 


tlh ] 
^41 


#32 


\wn Wi 


Wl 


Wl 


W2,W3, 


W5, W2, W4 


W5,W3,W2, Ws 


Ws 


11) I ] 


#33 


Iwn Wi 


U>2 


Ws 


W4 


W2 


W5 


Wo 


W2 


Ws 


Wo, 


Wl 


Ws 


Wo 


W4 


Wo] 


#34 


llDn 11) 1 


W2 


Wl 


Wl 


Ws 


W4 


Ws 


W2 


W4 


We, 


W4 


Wl 


W5 


W4 


11} q] 


#35 


\wn Wi 


W2 


W2 


Wl 


Wo 


W2 


W2 


Ws 


W4 


Ws, 


We 


W4 


Ws 


We 


We:] 




[Wo , W/J 


W2 


W2 


Ws 


Wo 


W4 


W4 


W4 


W2 


Ws, 


We 


W4 


W2 


We 


Ws] 


#37 


[lUO.Wl 


W2 


Ws 


Wl 


Wo 


Ws 


W2 


W4 


Ws 


We, 


W7 


Ws 


W4 


W7 


We] 


#38 


[wo.wi 


W2,W2, Ws 


W4 


W5,W5 


We 


W7 


Ws, Wg,We, Wj, Wg 


Ws] 


#39 


[wo.wi 


W2 


Ws 


Wl 


W4 


W5 


We 


W2 


Ws 


W7, 


Ws 


Ws 


Wo 


Ws 


Wg] 


#40 


[wo , Wi 


W2 


Ws 


W4 


Wo 


W5 


We 


We 


Ws 


Wj, 


Ws 


Ws 


W2 


Wg 


W7] 


#41 


[wo.wi 


W2 


Ws 


W4 


W5 


Wl 


We 


W7 


W4 


Wo, 


Ws 


W8 


We 


Ws 


Wg] 


#42 




W2 


Ws 


W4 


W5 


We 


W7 


Ws 


Wg 


Wio 


, Wii,Wi2, 


Wis 


W14 , Wl 
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21 



4 , 


Q 




r 

, 


1 


r 

, 


4 


r 

, 


Q 


r 

, 





r 

, 


4 


r 

, 
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Pattern number 



Generating function 
1 



1,2,3 



1 - s 



4,5,6,7 



8,9 



10,11,16,17 




12,13,14,15 



l + 3s 



{1-sf 
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18,22 
19,24 
20,23 
21,25 

26,30,31,33 

27, 29 

28,32 

34,36 

35 

37 

38 

39,40,41 

42 



{l-sf 
(1 + s) (l + 6s + 3s2_s3) 

(1 + s) (l+4.s-g^) 

l + 8s + 12s2 + i8s3 + 13s^ + 6s5-2s6 

(l + s)(l + s + s2)(l-sf 
l + lls + lls2-7s3 

l + 5s + 4s2_3s3_^g4 
(l + s + s2) (l + 7s-2s2) 

1 + 12 s + 16 + 3 s3 + 8 + 4 _ 2 

{l + s)il-sf 
1 + 12 s + 24 s2 + 34 + 27s^ + 14 
(l + s)(l + s + s2)(l-s)4 

l + 17s + 19s2-5s3 

l + 18s+lls2 + i0s3 + 8s^ + 8s5 
l + 23s+15s2 + 5s3_^4g4 

l + 41.', + 3.s- + 35.v' + 16s4 



(1 - sr 



Appendix C. Admissible patterns for q = 3 



#1 : [wo, 



#2 

#3 
#4 
#5 



[wo, 

[wo, 
[■Wo, 
[wo. 



#6 : [wo, 
#7: [wo, 
#8 : [wo, 
#9 : [wo. 



Wi,Wi,Wi, 
Wo, Wi,Wo, 

Wi,W2,W2, 
Wj , , , 

Wi,Wi,W2, 
Wi,W2,Wg, 
Wi,W2,Wi, 
Wi,W2,Ws, 



Wo,Wi, 
Wi,Wo, 
Wo, Wj, 

■W2, W3, 
W3,W4, 
W3,W4, 

Ws,W4, 

W4,W5, 



Wi,Wi,Wo] 
Wi,Wo, Wo] 

Wi,W2, Wo] 
W2,Wo,Wi] 
Wj, Ws, Wg] 
W2,W4,Ws] 
Wi,W5, Ws] 
W2,W4,W5] 
We, Wy, Ws] 
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Appendix D. Some algebraic invariants for g = 3 

The lower index is the pattern number. 

(a) _ (wi^W3 - Wo W2^) 



[1 

[6: 
[7: 

[9 

[10: 

[11 

[12: 
[is: 

[14 

[15: 
[16: 
[17: 
[is: 



W2 {Wi^ - Wq W2) {W2 Wi^ - Wq ' 



A5 



(6) _ (W2 Wl^ - Wowi^ 



Ws {Wi^ - Wq W3) (wi^Wg - Wq ^2^) " 



. (a) _ ^^1 ^2 

6 ~ 9 I 9 ■ 

W2 + Wl 



(6) _ {wiW2'W4 - WqWs'^)^ 



W3 {Wl W2 — Wo W3) {Wl W2 W3 — Wi'^Wo) ' 

{wi W2 W3 - Wi'^WoY 
W4 {Wi W2 — Wq W4) {wi W2 W4 — Wq W^"^) 



^ (g) _ Wq {W3 W5 - W4^) ^ (b) _ Wq (^3 W5 - W4^) ^ (c) _ ^2 W3 

{'W4W2 — WiWs) Wl [WiWi — W2W3) W2 Wl ^4 

References 

M.P. Bellon, J-M. Maillard, and C-M. Viallet, Infinite Discrete Symmetry Group for the Yang- 
Baxter Equations: Spin models. Phys. Lett. A 157 (1991), pp. 343-353. 

M.P. Bellon, J-M. Maillard, and C-M. Viallet, Infinite Discrete Symmetry Group for the Yang- 
Baxter Equations: Vertex Models. Phys. Lett. B 260 (1991), pp. 87-100. 

M.P. Bellon, J-M. Maillard, and C-M. Viallet, Rational Mappings, Arborescent Iterations, and 
the Symmetries of Integrability. Phys. Rev. Lett. 67 (1991), pp. 1373-1376. 

R.J. Baxter. Exactly solved models in statistical mechanics. London Acad. Press, (1981). 

F. Y. Wu, The Potts Model. Rev. Mod. Phys. 54 (1982), p. 235. 

J. Ashkin and E. Teller, Statistics of two-dimensional lattices with four components. Phys. Rev. 
64 (1943), p. 178. 

M. Kashiwara and T. Miwa, A class of elliptic solutions to the star-triangle relation. Nuclear 

Physics B 275([FS17[) (1986), pp. 121-134. 
S.Howes, L.P. Kadanoff, and M. den Nijs, Quantum model for commensurate-incommensurate 

transitions. Nucl. Phys. B 215([FS7[) (1983), p. 169. 
H. Au-Yang, B.M. Mc Coy, J.H.H. Perk, S. Tang, and M.L. Yan, Commuting transfer matrices 

in the chiral Potts models: solutions of the star-triangle equations with genus > 1. Phys. 

Lett. A123 (1987), p. 219. 
R.J. Baxter, J.H.H. Perk, and H. Au-Yang, New solutions of the star-triangle relations for the 

chiral Potts model. Phys. Lett. A128 (1988), p. 138. 
R.J. Baxter, On Zamolodchikov's Solution of the Tetrahedron Equations. Comm. Math. Phys. 

88 (1983), pp. 185-205. 

G. Falqui and C.-M. Viallet, Singularity, complexity, and quasi-integrability of rational 
mappings. Comm. Math. Phys. 154 (1993), pp. 111-125. 

M. Bellon and C-M. Viallet, Algebraic Entropy. Comm. Math. P hys. 204 (1000), pp. 42.5-437 



Gnu MultiPrecision. Arbitrary precision arithmetics software, see http://www.swox.com/gmp/. 
M.K. GizatuUin, Rational G-surface. Izv. Akad. Nauk SSSR. Ser. Mat. 44(1) (1080), pp. 110- 

144. English translation is in "Math. USSR, Izvestiya", vol. 16 (1981), no.l, p.103-134: MR. 
M.P. Bellon, Algebraic entropy of birational maps with invariant curves. Lett. Math. Phys. 50 

(1999), pp. 79-90. 

D. Hansel and J. M. Maillard, Symmetries of models with genus > 1. Phys. Lett. A 133 (1988), 
p. 11. 

M. van Hoeij. An algorithm for computing the weierstrass normal form. In International 
Symposium on Symbolic and Algebraic Computation (ISSAC '95), pages 90-95, (1995). 



A classification of four- state spin edge Potts models 



24 



[19] J-C. Angles d'Auriac, J-M. Maillard, and C-M. Viallet, Random Matrix Theory and higher 



genus intenrnhility: the quantu 



4801-4822. cond-mat/0205101 



n chiral Potts model. J. Phys. A: Math. Gen. 35 (2002), pp. 



[20] M.P. Bellon, J-M. Maillard, and C-M. Viallet, Integrable Coxeter Groups. Phys. Lett. A 159 
(1991), pp. 221-232. 

]21] S. Boukraa, J-M. Maillard, and G. RoUet, Integrable mappings and polynomial growth. Physica 

A 209 (1994), pp. 162-222. 
J22] H. Meyer, J.C. Angles d'Auriac, J.M. Maillard, and G. RoUet, Phase diagram of a six-state 

chiral Potts model. Physica A 209 (1994), p. 223. 
J23] H. Meyer and J-C. Angles d'Auriac, Phase diagram of a six-state chiral Potts model. Phys. 

Rev. B 51(10) (1995), p. 6379. 
J24] V.A. Fateev and A.B. Zamolodchikov, Self-dual solutions of the star-triangle relations in Zjv 

models. Phys. Lett. A92 (1982), p. 37. 
J25] M.P. Bellon, J-M. Maillard, G. RoUet, and C-M. Viallet, Deformations of dynamics associated 

to the chiral Potts model. Int. Journ. Mod. Phys. B 6 (1992), pp. 3575-3584. 



